Modeling f(R) gravity in terms of mass dilation rate 
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We review the conformal equivalence in describing the background expansion of the universe by 
f(R) gravity both in the Jordan frame and the Einstein frame. In the Jordan frame, we present the 
general analytic expression for f(R) models that have the same expansion history as the ACDM 
model. This analytic form can provide further insights on how cosmology can be used to test the 
f(R) gravity at the largest scales. Moreover we present a systematic and self-consistent way to 
construct the viable f(R) model in Jordan frame using the mass dilation rate function from the 
Einstein frame through the conformal transformation. In addition, we extend our study to the 
linear perturbation theories and we further exhibit the equivalence of the f(R) gravity presented in 
the Jordan frame and Einstein frame in the perturbed space-time. We argue that this equivalence 
has solid physics root. 
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I. INTRODUCTION 



There has been conclusive evidence indicating that our 
universe is experiencing an accelerated expansion. This 
acceleration is believed driven by a so called dark energy 
(DE) in the framework of Einstein's general relativity. 
The simplest explanation of such DE is the cosmological 
constant. But it suffers serious problems such that its 
value is far below the prediction of any sensible quantum 
field theories and it inevitably sustains the coincidence 
problem with the mysterious same order energy density 
as the matter field today. 

There exists an alternative way to explain the accel- 
eration of the universe expansion by modifying the Ein- 
stein's gravity. One simplest attempt is called f(R) grav- 
ity, in which the scalar curvature in the Lagrange density 
of Einstein's gravity is replaced by an arbitrary function 
of R [1]. However it is quite non-trivial to construct 
the viable f(R) model which can satisfy both cosmo- 
logical and local gravity constraints |2h10||. Only a few 
viable / (R) models with simple analytic forms have been 
found flfj| - tl3| . However, if we do not limit to these f(R) 
models with explicit expressions, in a more general case, 
f(R) gravity in cosmology can leave enough freedoms to 
accommodate any desired expansion histories of the uni- 
verse [14j [15[ . A general phenomenological f(R) cosmo- 
logical model were formulated numerically by specifying 
the effective DE equation of state w and the present value 
of the parameter B = j^§ [H QJ]. 

On the other hand, it is always possible to transform 
the theory of f(R) gravity presented in the Jordan frame 
to the Einstein frame by using the conformal transfor- 
mation [17] . In the Einstein frame, the model contains 
the coupling between the matter and an additional scalar 
field [18[ . If the scalar field plays a role of DE in the uni- 
verse, DE would interact with DM. The phenomenologi- 
cal DE and DM interacting model has been studied ex- 
tensively in literatures, see [19l425| and references therein. 
More interestingly, it was found that a general f(R) grav- 
ity can be consistently constructed in a covariant form in 



terms of the mass dilation rate function in the Einstein 
frame[2g. This mass dilation rate function marks the 
coupling strength between DE and DM. The condition 
for the f(R) cosmology to avoid the instability at high- 
curvature region and be consistent with CMB observation 
require the energy flow from DE to DM in the framework 
of the interacting model [26[ . This was found helpful to 
alleviate the coincidence problem J19l |23| . 

In this paper, we will further examine the equivalence 
of the f(R) gravity presented in the Jordan frame and 
the Einstein frame to describe the cosmology. In the 
background Jordan frame, we will derive the analytic 
f(R) functions to give the same expansion history of 
the universe as the ACDM model. This can provide fur- 
ther insights on how cosmology can be used to test the 
f(R) gravity at the largest scales. In addition to model- 
ing f(R) gravity in Jordan frame using the effective DE 
equation of state, we will present a systematic and self- 
consistent way to construct the f(R) model in terms of 
the mass dilation rate function in Einstein frame. More- 
over we will develop linear perturbation theories to fur- 
ther exhibit the equivalence of the f(R) gravity in the 
Jordan frame and the Einstein frame. 

This paper is orgnized as follows: In section [Hi we will 
discuss the f(R) gravity in the background cosmology. In 
the Jordan frame, we will derive the general analytical so- 
lutions for the f(R) model which has the same expansion 
history as the ACDM cosmology. In the Einstein frame, 
we will argue that once the mass dilation is specified, 
the f(R) gravity can be constructed accordingly. In sec- 
tion IIII1 we will extend our discussion to the perturbed 
space-time. We will show that perturbation theories in 
Jordan frame and Einstein frame can be consistently con- 
nected by conformal transformations. In section HVl we 
will construct a viable f(R) cosmological model in terms 
of mass dilation rate specified in the Einstein frame and 
discuss the subhorizon approximation of the constructed 
f{R) model. Our conclusion and discussion will be pre- 
sented in the last section. 



II. BACKGROUND DYNAMICS IN f(R) 
COSMOLOGY 



In the Jordan frame, the matter field is conserved, which 
satisfies 



We start with the 4-dimensional action in f(R) grav- 
ity 



S = Y2 I d A x^f{R) + I d 4 xC (m) 



(1) 



where k 2 = 8ttG. The field equation can be obtained by 
varying the above action with respect to g^, which leads 
to 

FR va ,-\fg^-V^ v F + g^UF = ^T^ , (2) 

where F — ^4 and 



DF = 



-MV^gg^d^F) 



(3) 



If we define the left hand side of Eq. @ as a tensor E^ 



^ = FR llv --fg ia ,-'S7 ll 'S7 v F + g lu ,DF , (4) 

we can recast the modified Einstein's equation into a sim- 
ilar form to the standard Einstein equation 

E^ = « 2 T(") . (5) 

The trace of the above equation gives rise to a scalar 
equation 



1 1 K 2 , 

f = '±UF+-FR-—T^ 
J 2 2 2 



(6) 



Taking the derivative of / and further noting that 

V M / = FV M i? , (7) 

we obtain 

V Al i?=-^V M T( m ) + |v M DF + i?V A1 ln J F . (8) 

We can see that F is a free scalar, which can be used to 
characterize the f(R) gravity. 

A. The Jordan frame 

In the Jordan frame, we consider the universe 
which is described by the fiat Friedmann-Robertson- 
Walker(FRW) metric 



ds = —dt + a dx 



2 , „2 7 2 



(9) 



From Eq.([5]), we can obtain the equation governing the 
dynamics in the background spacetime 



H 



, FR- f TT F : k 2 



liF H F + 3F Pn 



(10) 



3Hp m = , 



(11) 



where dot denotes the derivative with respect to the cos- 
mic time t. 

From Eq. (fT0|) . / can be presented as 



/ = FR- 6FH 2 - 6iJ 2 -^ + 2k 2 p„ 
dx 



(12) 



where x = hi a. 

Taking the derivative of the above equation and further 

df_ _ pdR 

dx dx ' 



noting that £ = F4§, we obtain 



d 2 „ ,ld\aE ^dF d\nE 3n°,e" 

-F4 



3.i' 



-F + ( 
dx 2 K 2 dx 



where 



dx dx 



E 



, 
(13) 



E = 



IP 



9. 



K 2 n° 

_ K Pm 



3H 2 ' 
,dE 



(14) 
(15) 



R is in the unit of Hq and k 2 = 1 . 

For illustrative purposes, we take an expansion history 
in the Jordan frame that matches a DE model with equa- 
tion of state w, 



e = n° m e- 3x + (1 - nl l )e' i ^ {1+w)dx 



(16) 



In [15| a family of numerical solutions of f{R) models 
that match the representative expansion histories of DE 
model was obtained. Here we find that for the f(R) 
model to mimic the same expansion history as that of 
the ACDM cosmology in the background E = 9,^ n e~ 3x + 
(1 — Q^J , we can analytically solve Eq. (fT3| and obtain 
the general analytical solution 

F(x) - 1 = 
C{e^r-(l-^r-(Sf m )- p - 2 F 1 [q. ) p.,r.,z] 

+D(e Sx r+(l~n° n r+(n m )- p+ 2 F 1 { q+ ,p + ,r + ,z] , 

(17) 

where 2-F\(a, 6; c; z) is the hypergeometric function, 
which is defined as 



2 F 1 (a,b;c;z) = ^ 



n=0 



(fl)n(b)n Z n 

(c)„ n! 



(a) n = a(a + 1) • • • (a + n — 1) 
C and D are coefficients which will be determined by 



boundary conditions. The indexes in the solution read, 
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p+ - 


12 




5- V73 


P- - 


12 


q+ = 


1 + V73 
12 


q- = 


I-a/73 

12 


r+ = 


= 1 + f , 

6 


r_ = 


_ 1 \/73 



(18) 



and 



_3x 



(i 



,,3.r 



30i 



i?- 12(1-00.) 



(19) 



Since the index p_ is negative, the fist term of Eq.pTj) 
will diverge when x goes to the minus infinity, which is 
not acceptable in physics. Besides, a viable f(R) model 
should be in a form of "chameleon" type [28j, which 
means that in the high curvature region it should go back 
to the standard Einstein gravity to pass the local Solar- 
System tests. The boundary condition for the scalar field 
F(x) thus requires 

lim F(x) = 1 , (20) 

x— > — oo 

and we need to set C = 0. The solution turns out to be 
F(x) = 1+D(e 3 *r+ (l-0 P+ (0^ + 2^ib+,P + ,r + ,z] 



(21) 
D above is a free parameter which characterizes the dif- 
ferent f(R) models which have the same expansion his- 
tory as that of the ACDM model. 

With the solution of F(x), we can easily figure out the 
analytical expression for f(R) by doing the integration 



f(R) = f(R(z)) = J F(z)™dz . 

We obtain 

f(R) = [L A - 6(1 - 0° n )](l - e) + Constant 
where 

e = D^p\(l-a° m r+(a m )-^(e 3x r+ 2 F 1 [ q+ ,s + ,r + ,z] 

r E[P+\ 



(22) 
(23) 



and 



La = R-2A = R-6(l-n° n ) 

-7+\/73 



f(R) 
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FIG. 1: f(R) models that have the same expansion history as 
the ACDM cosmology. 



r# is the complete Euler Gamma function, which is de- 
fined as 



r*(*) = r 

Jo 



(25) 



t^'e-'dt 



When D = 0, Eq.(|23|l should be consistent with Eq.(|T2|) 
and should go back to the ACDM model. This puts the 
constant in Eq. (|2"3"]) into 6(1 — O^J. Thus finally we have 
the analytic f(R) form 

/(A) = [£ A - 6(1 -<&)](! -e) + 6(1 -(&) . (26) 



s+ 



12 



Eq. ([2^|) is a general analytical expression for f(R) 
models that can have the same background expansion 
history as that of the ACDM cosmology. The solution 
has only two free dimensionless parameters, namely, 51^ 
and D. In order to avoid the short-timescale instability 
at high curvature [15|, it requires D < 0. Moreover for 
a viable f(R) model, it also requires F > 0, fuR > 
[18| . In Fig [TJ we illustrate the f(R) models obtained 
from Eq. ([26"]) and the blues lines indicate the viable f(R) 
models which satisfy the condition mentioned above. 

If we do not start from the concrete model describ- 
ing the expansion history of the universe, there is a lot 
of freedom for choosing the scalar field F(x) except the 
boundary condition Eq. (|2"0)) . The f(R) gravity can be 
constructed only when the scalar field F(x) can be rea- 
, sonably chosen based on physical motivations. However, 
in the Jordan frame, there is no clear physical motivation 
to provide a particular form for the field F(x). We need 
to turn to the Einstein frame. In the Einstein frame, 
F(x) relates to a scalar T(x) represen ting the mass di- 
lation rate d ^y" in gravitational field [lj, [26[ and T(x) 
(24) can be deemed as the coupling strength in the framework 



of interacting dark energy models [H IM IM IH • F{x) 
can be determined by constructing Y which, in turn, the 
f(R) model can be constructed. In the next section we 
will shift our attention to the Einstein frame. 



B. The Einstein frame 

The equations of motion in the Einstein frame can be 
obtained by using conformal transformations 



H 



1 , ldlnF x 



(27) 



As shown in [26] , after the conformal transformation the 
expansion history in the Einstein frame is governed by 



dx QH 




(28) 


dR dY 9 
-rrr + 36r— + 72Y 2 = 
dx dx 


p 

-3K 2 p m (-~- + 1) 
H 


. (29) 


dpm „ - r _ 

dx +6pm ~ jj Pm ' 




(30) 



where Y is the coupling strength between the scalar field 
and matter which can also be explained as the mass di- 
lation rate function [2a|. From the relation between F 

l~2 



and r and expressing F into <p, InF = K\/^tp, we can 
rewrite the continuity equation for the scalar field into 



dip 
dx 



Vqy 

kH ' 



(31) 



Once the dilation function is specified, we can explore the 
expansion history of the universe in the Einstein frame. 

r has clear physical meaning, which represents the cou- 
pling strength of the scalar field to the matter field. Al- 
though T is a free field, the choice of Y must be consistent 
with the viability condition of the f(R) gravity. A viable 
f(R) model must pass the stringent local test, which re- 
quires that when R —> oo, the model should go back to 
the standard Einstein gravity, so that limj^._ 00 Y = 0. 

Eqs. (f28|) - (l30|) governing the expansion history are not 
stable when T < at high curvature region, except that 
the coupling strength is extremely small, for example 
|r| < 10~ 10 . This phenomena is consistent with that 
discussed in [29]. Thus a viable f(R) model requires 
that r is positive and drops to sufficiently small values 
in the early time of the Universe. 

r also is the dilation rate of the mass of test particles in 
gravitational fields. To satisfy the weak equivalent prin- 
ciple, r should be independent of the species of matter 
[301 ] . One natural choice of the Y form is to consider it as 
a geometrical quantity, a function of R. There is quite 
a wide range of the choices for Y which can satisfy the 
above conditions. In this work, we only present one of 
the simplest choices, where we describe T as a function 



of the scalar Ricci curvature R in the Einstein frame as 

a 



Y = -^- 



R + /3 



(32) 



where a and /3 are constants. Eq. (|32p is just one kind of 
possible phenomenological descriptions since we are lack 
of the knowledge on the nature of the dilation rate. 

Once we specify the dilation rate Y, we can get the 
f(R) gravity in the Jordan frame by using the conformal 
transformation from the Einstein frame to the Jordan 
frame. 



III. PERTURBATION THEORY IN f(R) 
COSMOLOGY 

In this section, we will extend our discussion to the per- 
turbed space-time. We will first review the cosmological 
perturbation theory in f(R) gravity in the Jordan frame. 
We will show that perturbations can be consistently con- 
nected between the Jordan frame and the Einstein frame. 



A. The Jordan frame 

We focus on the scalar perturbation. The perturbed 
line element in Fourier space reads, 



ds 2 



M, 



= a 2 [-{\ + 2iPY (s >)dT 2 + 2BY i {s 'drdx l 

+ (l+2<f>Y^)6 ij dx i dx j +EY l { J s) dx i dx 3 ] ,(33) 



where tp,B,(p,E are scalar fields which indicate the per- 
turbations. Y^ is the scalar harmonic function. The 
perturbed form of the modified Einstein equations Eq.([5]) 



can be obtained by writing 

where 

SZp, = dFR^+FSR^-^Fg^SR- 



(34) 



/#9W 



- SV^VuF + Sg^DF + g^SDF , 

and expressing 

ST (m)0 °Y^ = -S Pm Y^ , 
5T (m)i °Y} 8) = (p m + Pm )(v m + B)Y} s) 
ST {m]0 l Y^ = -( Pm + Pm )v m Y^ i , 



(35) 



ST, 



(m)j 



iy( s )'' _ 



(Sp m S) 



n m /)^. 



(s)i 



(36) 



Inserting the line element Eq. (|33[) into the above equa- 
tions, we can get the perturbation form of the modified 
Einstein equation. 



From (5E = K 2 STf m \ we obtain 



K 2 1 

— 5 Pm a 2 = F[-k 2 cf) + in{'Hil)~cj) l ) + knB~-k 2 E} 
+ F'(^B + m^-\<f>')-\nHF 



3 n" 3 k 2 

2 a 2 2 



(37) 



where the prime denotes the derivative with respect to 
the conformal time r. 

The trace part <5IV = K 2 8T^ m ^ gives 

K 2 a 2 S Pm = F\--k 2 ^--k 2 <j)--k 2 E + A—4, + -kB' 
3 3 9 a 3 

+ \kUB - 2H 2 ip + 2Htp' - AHtf - 20"] 
+ F'[^kB + i;' + 2nip-2<f)'] 



After inserting Eq. (grj into Eqs. (|57 ]) -(|3n ]) . we find that 
Eqs. (|37|) -(|40 |) keep the same forms under the infinitesi- 
mal transformation, which shows that Eqs. (f3"T |) - f4T))) are 
covariant. Furthermore, if we take F — > 1,SF — > 0, Eqs. 
(I57l) - P0)) can go back to the standard cosmological per- 
turbation equations in the Einstein gravity as presented 
in P. 

After we get the perturbed form of the modified Ein- 
stein equations, we turn to derive the equation of motion 
for matter field in the Jordan frame. From 



TTJVrp{m) 

v ■*- fJ,V 







(42) 



we can obtain its perturbed form SVT^™' = 0. Looking 
at the zero-th component, we have 



5' m + m{ 



Sp n 

Sp„ 



W m )5r, 



6F[H 2 



a 
a 



-k 2 } - 5F" - 5F"H + 2</>F" 



-(1 + w m )kv m - 3(1 + w m )4> 



(43) 



(38) 



The i-th component gives 



and the anisotropic part 5Tii J = k STi m -\i J (i ^ j) gives 



K a Pmttn 



F[-k 2 %l) - k 2 (, 



2kUB + nE' - -k 2 E 
6 



-E" + kB'} + F'{-E' + kB) - k 2 SF. 



(39) 



The last part ST, l = K 2 <5T( m ) l leads to 

1 



K a (p m +p m )Vr, 



F[k(/>' - kUip + 2H 2 B B 

+ -kE'} +F'[HB - -kip] 

+ -kSF' - -kHSF - -BF". 
2 2 2 

(40) 

The above perturbation equations are covariant. This 
can be easily seen by investigating the infinitesimal coor- 
dinate transformation, 

aff = x ^ + Sx^ , 
Sx° = C°Y is) , 
Sx* = (3Y^ . 

We use the symbol with "hat" to indicate the quantities 
after coordinate transformation. Under the infinitesimal 
transformation, we can show that the perturbed quanti- 
ties in the line element Eq. (|33|) behave as 25 1 



[(pm +Pm)(B + V m )} + 4H(p m + Pm)(B + V m ) 

2 
= [Pm + p m )kip + kSp m - -kp m IL m . (44) 

Similarly, it can be shown that the above equations are 
covariant under the infinitesimal coordinate transforma- 
tion. 

The freedom of choosing coordinates in general rela- 
tivity will lead to ambiguity in defining a real density 
perturbation from the coordinate transformation. The 
systematic way to get rid of the ambiguity is to fix the 
gauge, which has been well discussed in [3lJ. Here we 
will argue that the gauge conditions in the Jordan frame 
are not always the same as the conditions in the Einstein 
frame after the conformal transformation. 

We only focus our discussion on the most commonly 
used gauges, the Newtonian gauge and the Synchronous 
gauge. The Newtonian gauge is defined by setting B = 
E = and these conditions completely fix the gauge 



B 


E' 
'2k 2 


E 




2k 





> 7- rw9 



/3 



(45) 



The perturbations in this gauge can be shown as gauge 
invariant. For instance, the matter density perturbation 
can be presented as 



^F« = (ij) - ?' - U£?)Y^ , 
4>Y& = 



5 (N) = S-^ 
P 



'<! 



2k 2 ' 



(46) 



BY, 
EY 



W = 
M _ 



(4>-hp-ue)Y [s) 

(B - k£° ~ (3')Y} S) , 



(E + 2k/3)Y t 



(«) 



(41) 



where 5 is a general perturbation. We can show that S^ 
is invariant under the coordinate transformation Eq.(l4Tj). 
However, in contrast to the Newtonian gauge, the Syn- 
chronous gauge is not completely fixed because the gauge 



condition ip = B = only confines the gauge up to two 
arbitrary constants C± , C2 , 



e - 


9i 

1 
a 




- 


■ - kc \i 


f dr 

a 



Co 



(47) 



In order to fix the Synchronous gauge, additional con- 
ditions are called for. Usually, Ci is fixed by specifying 
the initial condition for the curvature perturbation in the 
early time of the universe and C\ can be fixed by setting 
the peculiar velocity of DM to be zero, v m = 0, which is 
equivalent to say that the gauge is at rest with respect 
to the cold DM and the observer is comoving with cold 
DM. After fixing C\ , C2, the Synchronous gauge can be 
completely fixed. 

As long as the Synchronous gauge is completely fixed, 
it will be equivalent to the Newtonian gauge. For in- 
stances, the density perturbations in different gauges are 
related by 



§(N) =S (S) + Prn'H 



N y S 



k 



(48) 



F is a free scalar field which characterizes the deviation 
from the Einstein gravity. The perturbation 5F can be 
obtained by 



dF 

SF =m SR 



(53) 



where 



a 2 SR = H(18& - m>') + 4fc 2 $ + 2kH - 12*— + 6$ . 

(54) 
Using Eq. ([54]) and Eq. (|4"9")l to eliminate $ and fc 2 <£ in 
Eq.([5"0"|). we obtain 



k 2 
SF" + 2USF' + a 2 (— + M 2 )5F 



-a 2 Sp m + F'{AU^ - 3$' + *') + 2^F" ,(55) 



where 



M 2 = \(J--R) , 

o JRR 



(56) 



where v N = v — §b, v s = 0. 

We will show later that the gauge condition for the 
Newtonian gauge will be the same both in the Jordan 
frame and the Einstein frame. However we will see 
that the gauge condition for the Synchronous gauge will 
change when we transform the Jordan frame to the Ein- 
stein frame by using the conformal transformation. 

Since the Newtonian gauge has this advantage, we will 
work in the Newtonian gauge hereafter. We use the 
Bardeen potentials [32J $ = <f>, * = ip to represent the 
space time perturbations. Since we are only interested in 
the DM dominated period in the f(R) cosmology, we set 
p m = and Sp m — 0. The perturbed Einstein equations 
can be reduced to 



and 



K s 2 

^-op m a 


= 


F[-fc 2 $ + 3-H(-H*-$')] 






+ 


F'(m^ - -$') - -n 2 SF 






+ 


3 n" 3 k 2 

1^SF--HSF'~-SF 
2 a 2 2 


(49) 


SF" 


= 


2 2 a" 
-H6F' + F[ k 2 ^ - -fc 2 $ + 4—* 

3 3 a 




- 


2"H 2 * + 2W$>' - m& - 2$' 


"'] + 2*F" 




+ 


F'[*' + 2-H*-2$'] 






+ 


»!«■+£- §*1 , 


(50) 


SF 


= 


-*-$ , 


(51) 


1 2 2 

k a p m v m 


= 


F[k& - kW&] - -kF'y 





+ -kSF' - -kUSF 



(52) 



D _ 6a" _ dF 

H = ~^T 1 JRR — TT5 
a^ oR 



(57) 



Eq. (l55")) is consistent with Eq.(8.92) in [Ig, if we express 
it in the cosmic time t. 

^From the equations of motion Eqs. (|4"31 (|44p , we obtain 



S' m = -kv m - 3$' 

vL + Uv m = fc* 



(58) 
(59) 



Combining Eqs. (1581 159]) . we have the second order differ- 
ential equation for the evolution of the matter perturba- 
tion 



US' + fc 2 * = -3H& - W 



The Einstein frame 



(60) 



We will first derive the perturbed form of the equation 
of motion for matter fields Eq. (14"2"|) under the general con- 
formal transformation. Then we will fix the freedom of 
the general conformal transformation by setting tt 2 = F 
to obtain the perturbation form of the modified Einstein 
equations. We will show that the equations obtained un- 
der conformal transformation from the Jordan frame are 
consistent with the results derived by doing perturbations 
directly in the Einstein frame. 

In the background, the conformal transformation is de- 
fined by re-scaling the line element of the Jordan frame 



Is 2 = n 2 ds 2 



(61) 



However, in perturbed spacetime, we need to take ac- 
count of the perturbation of n. And thus the full line 
element reads 

ds 2 = (n + SQ) 2 ds 2 

= & 2 [-(l + 2^)rfr 2 +2Bdrdx i 

+ (l + 24>)S ij dx l dx j +Edx i dx j ] . (62) 

Since we are only interested in the linear perturbation, 
we can expand the above expression to the linear order 



i> -- 


-- i>- 


-tflnfi 


<$> -- 


-- 4>- 


■SlnQ, 


B = 


= B 


1 


E = 


= E 





(63) 

where the symbols with "tilde" indicate the quantities in 
the Einstein frame. 

For the clear discussion, hereafter we will use "tilde" to 
indicate quantities that are different in the Einstein frame 
from their values in the Jordan frame. For the quantities 
without "tilde", they are the same in both frames. We 
will use "hat" to indicate the coordinate transformation 
which should not be confused with the conformal trans- 
formation. The differential operator d does not depend 
on the geometric structure g a b on the spacetime mani- 
fold, it will not change under the conformal transforma- 
tion so that the ordinary derivative operator d a and the 
conformal coordinates dr, dx will be the same in both the 
Einstein frame and the Jordan frame. 

Noting 

<nnQ = <51nf2-(lnfi)'£° , (64) 

and considering Eqs. (|63l41l) . we can show that, under the 
infinitesimal coordinate transformation, the perturbed 
quantities ip, cj) behave as 



where 






H = H- (Inn)' 



(65) 



(66) 



ft is clear that in the coordinate transformation ip, <fi 
behave in a similar way as in the Jordan frame. This 
point also holds for other perturbation quantities 

B = B-k^-0' , 

E = E + 2k/3 . 

Now we can discuss the gauge conditions in the Ein- 
stein frame. In the Newtonian gauge, the gauge condi- 
tions B = E = fix the gauge as that in the Jordan 
frame, namely 



(67) 





B 


E' 




~T~ 


~2k~z 


p = 


E 
2k 





Thus the Newtonian gauge conditions keep the same un- 
der the conformal transformation. 

However, in the Synchronous gauge, after the confor- 
mal transformation, the gauge conditions in the Einstein 
frame read 



ip = 6 In , 
B = , 

V m = . 



(68) 



We will find that v m = is a physical choice because it 
is the solution of the equation of motion in Eq. (|75f . It is 
clear that in the Einstein frame the gauge conditions for 
the Synchronous gauge are not exactly the same as those 
in the Jordan frame. 

Considering the gauge conditions for the Newtonian 
gauge are the same in both the Einstein frame and the 
Jordan frame, we will work in Newtonian gauge hereafter. 

In the background, the conformal transformation leads 
the equation of motion of the matter field Eq. (|42l) in the 
Einstein frame to the form 



where f ^ = ~g^ T ~g va T Tai T^ = g^ T g ua T T<7 and 

1 



rp ~*- rp rp(j,U 



I T 

rpfXU rp 

n 6 ' ~ ft 4 ' 



(69) 



(70) 



_rpd_n - 1S ^ e C0U pi m g vector. We only consider the 
evolving universe where the coupling vector is time-like. 
We focus on the DM dominated phase with p m = 0. 
The conservation equation Eq. (jM)) can be recasted into 



(71) 



OT 

where U° = 1/a, T is the mass dilation rate 

r = - r0 91nn = 1 ppdm 
dr ifi dr 

The perturbed equation of Eq.([7Tj) reads 

s' m = -kv m - 3$' + ar* + a<5r , (73) 

where we have used 



(72) 



a a 

5T = M°(\nny -u°(s\nny 



(74) 



and SU° = -ipu . 

In the Einstein frame, the perturbation equation for 
the peculiar velocity has the form 



v' m + Hv m = aT(v t - v m ) + k& , 



where 



v t = k 



51nn 

(Inn)' 



(75) 



(76) 



Egs. ([751 175)) are consistent with the results in [23[ and [20| where we have used 
by setting the general coupling vector as Q^ = Tp m U^ , 



where U£ is along the direction of d v Vt 

Egs. (1731 [75]) can also be consistently obtained from 
the perturbation equations of Egs. (|551 15§f by doing the 
conformal transformation 

5 

° = ? l ' 

u° = nu° , 

pm ^ Pm i 

(5pm = (50 4 p m + Vt A 8p m , (77) 

<5 m = <5 m + 4(<51nQ) , 
* = *-51nrj , 
$ = $ - 5 In ft . 

Similarly, combining Eqs. ([73"l ITS")) , we can get the sec- 
ond order differential equation 

S' m + {H + aT)S' m + k 2 ^ 
= (H + aT)aST - kaTi t + (aST)' 
+ {H + aT){aT^ - 30') - 30" + (ar^)' . (78) 

Next, we consider the perturbed Einstein equations in 
the Einstein frame. We set ft 2 = F to fix the freedom 
of the general conformal transformations. Inserting the 
conformal transformations Eg. (1771) into Eqs. (|4"9"|) - (|52")) , we 
obtain 



-^-{5pm + 8pd)a 

2-2;- 

K a dp d 



-k 2 § + m{U^ -$') , (79) 

--fc 2 *- -fc 2 $ + 4— * 
3 3 a 



- 2H 2 * + 27i#' - 4H$' - 2$" 

o = -*-i , 

fc$' - kH^ = --n 2 a 2 [p m v m + (pd+Pd)vd\ 

where Vd, Pd and pd are defined as 

SlnF 



I'd 



-k- 



2aT 

3r 2 ^ 

Pd = — + v 

sr 2 

FR - f 



V = 



2k 2 F 2 



(80) 



The perturbed Einstein equations give rise to the pertur- 
bation of SV as 



5T = -H>T- — (SlnF)' . 
2a 

If we directly perturb this V expressed in Eq. (150)) 
kHv _1SFR (FR-f)SF 



2 F 2 



F 3 



after inserting 



R = 6F(3*r + T 2 + ^ + -r) 
a a A a 



(82) 



(83) 



(84) 



eliminating / by the Friedmann equation and further em- 
ploying 



14 l n" 

2 ~ a n o W CT^2 

k p m = 4-r^ - 2— - 6r 



(85) 



we can get the result in consistent with Eq. (|5Tj) . 

SF can be determined by Eq. ([55|) . Inserting Eq.([77|). 
(f55|) turns out to be 



SF" = \-k 2 + -K 2 a 2 p m -2a 2 T 2 + 2 — 
3 a 



1 a 2 



1 



45T' + 12ST - --^— ]<JF + ~K,a 2 FSp m 

3 JRR 3 

- (2"H + 4:aT)8F' - 2aFT#' - 4F51T' 

- 12F*aT + 6Fa$T , (86) 

where f RR is given by 



_ F' _ 2&FT 

fRR -ft- -gT ' 



(87) 



and i? can be obtained from Eg. (1541) in the Einstein 
frame. 

After fixing SF, we can obtain the consistent covari- 
ant perturbation for T by Eq. ([82|) . since the infinitesimal 
coordinate transformation for ST presented in Eq.([82|) 
satisfies the rules of the covariant transformation for a 
scalar field 



I A) 



ST = SV - Y'£ 



Similarly the perturbation for H , SH is given by its 
covariant configuration 



and so does ST 



3H = V^[/ M 



r = -u^Vfiinii 



(89) 
(90) 



IV. SUBHORIZION APPROXIMATION 



K 2 ti 2 8V = (SlnF)(9a'T + 2U 2 - 35 2 r 2 - — + 3ar' 



From the above analysis, f(R) model in the Jordan 
(81) frame is conformally equivalent to the interaction model 



in the Einstein frame. This equivalence holds not only 
in the background dynamics but also in the perturbed 
spacetime. Therefore, we can construct f(R) models by 
choosing a peculiar form of the mass dilation rate func- 
tion r in the Einstein frame. One of the simplest choice 
of the r is 

r = -^ 



1.6 



R + /3 

After fixing the mass dilation rate function, we can figure 
out the evolution of the cosmological model in the Ein- 
stein frame and then we can transform it into the Jor- 
dan frame by using the conformal transformation. In the 
lower panel of Fig [51 we present some examples for the ef- 
fective DE equation of state w for the f(R) cosmology in 
the Jordan frame constructed by selecting the above mass 
dilation rate function in the Einstein frame after the con- 
formal transformation. We can see that the effective DE 
equation of state in the constructed f(R) model presents 
the phenomena in the vicinity of —1 and with the —1 
crossing behavior, which is consistent with present cos- 
mological observations. This shows that with the appro- 
priate choice of the mass dilation function, which refers to 
the mass, a more fundamental quantity, we can construct 
the viable f{R) model to produce the observationally al- 
lowed effective DE equation of state. 

We have also shown that the linear perturbation equa- 
tions are exactly equivalent at all scales in both the Ein- 
stein frame and the Jordan frame. In the usual f(R) 
model, it was argued that the scale factor described by 
the f(R) cosmology during the matter phase behaves in- 
consistently with cosmological observations 9|. It is of 
interest to examine the f(R) model constructed from 
the peculiar mass dilation rate function we have chosen 
and see whether our f(R) model is viable if comparing 
with observations in the matter dominated phase. We as- 
sume that the quantities defined in the Jordan frame are 
consistent with the actual observational quantities. We 
will study the subhorizion approximation in the Jordan 
frame. 

Following [3, from Eqs. (|49l5ip . we obtain 



$ 
* 



-2 6lnF+ ¥2F 6p " 



1 



S\nF 



k 2 2F 



8p ri 



where SF is governed by Eq.(|55p 
S\nF k 2 1 



8p„ 



3F 



k 2 



■M 2 



(91) 



(92) 



and M 2 is defined by Eq. (|56|) . From Eq. ([60|) . the equa- 
tion for the growth function G(z 7 k) — d l ^ x D in the sub- 
horizon approximation reads 

^ + G 2 + (2+^)G-~O m/ ,(M) = , (93) 



i i i i I i i i i I i i i i I i i i i I i 



T 




k=0.01hMpc" , ,r o =0.2,a=1 
k=0.05hMpc" 1 ,r o =0.2,a=t 
k=0.1hMpc" 1 , r o =0.2,a=1j 
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FIG. 2: The evolutions of /i and the effective DE equation 
of state w in the Jordan frame for the f(R) models we con- 
structed. 



where D(z, k) 



8 m (z,k) 
S m ( Zi ,k) 



and 



H(k,x) = 



1 A + 3M 2 a(x) 2 /k 2 
F3(l + M 2 a(x) 2 /k 2 



(94) 



The M 2 can be very large in the early time of the uni- 
verse, since T —> as R — !• oo. When M 2 > p-, the 
regime is called "General Relativistic(GR) regime" where 
the matter density perturbation has the standard evolu- 
tion since the modified gravity doesn't show up in this 
regime where SF — > and $/^ — > — 1, /i(fc, x) — > 1. How- 
ever, in the late time of the Universe, when M 2 <C fr, 
the modified gravity effect has prominent influence on the 
structure formation because the ratio of $/"$> — > —1/2, 
which has been altered greatly from the standard Ein- 
stein model and in this case /i(k, x) — > -4- where Fq is 
the value of F today. This regime is called the "scalar- 
tcnsor(ST) regime" and obviously matter perturbation 
does not have the standard evolution as in the Einstein 
gravity in this regime. As shown in the upper panel of 
Fig[21 in the transition phase of these two regimes, /x(fc, x) 
is strongly scale dependent on the wave number of k. The 
transition of these two regimes can be characterized by 
M 2 = k 2 /a 2 which can occur as early as in the matter 
dominated epoch, even in small scales, in certain types of 
f(R) models in the expansion history of the Universe [2]. 
In our model, since the dilation rate function T will drop 
rapidly in the early time of the Universe, the transition 
usually happens very late z < 30 for interested wave 
number k . Thus in our constructed f(R) models, we 
have the standard matter dominated phase and matter 
power spectrum, which can avoid the problems for the 
usual f(R) model mentioned in Q. 
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V. CONCLUSIONS 

In this work, we have further disclosed the equivalence 
of the f(R) gravity presented in the Jordan frame and 
in the Einstein frame from the background dynamics to 
the perturbation theory. In the Jordan frame we have 
derived the general analytic solution for the f(R) model 
to have the same expansion history as the ACDM cos- 
mology. This analytic solution of the f(R) model can 
provide further insights on how cosmology can be used 
to test the gravity at the largest scales. Moreover we 
have presented a systematic and self-consistent way to 
construct the f(R) model in the Jordan frame through 
conformal transformation by using the mass dilation rate 
function in the Einstein frame. We have also developed 
linear perturbation theories to further exhibit the equiv- 
alence of the f(R) gravity in the Jordan frame and in 
the Einstein frame. We have shown that the f(R) model 
constructed from the mass dilation rate can produce ob- 
servational allowed effective DE equation of state and 
provide the matter phase in consistent with cosmological 
observations. It is of great interest to confront our f(R) 
model with latest observations in the future study. 

It is interesting to note that our f(R) model was con- 
structed from the mass dilation rate, which refers to a 
fundamental quantity, the mass, rather than the myste- 
rious DE equation of state. This gives the root for under- 
standing the equivalence between the f(R) model in the 
Jordan frame and the Einstein frame, since we know well 
that the concept of mass plays a central role in gravi- 
tational theories [13, l33|- The mass can be classified 
into the inertial mass and the gravitational mass. In the 
Jordan frame the f(R) gravity can be effectively treated 
with the modified gravitational constant G e // = % if we 
compare it with the Einstein gravity. This is effectively 
equivalent to rescale the gravitational mass rh g = my/F. 
However the inertial mass in the Jordan frame is con- 



served so that the equation of motion for a free particle 
in the Jordan frame is described by u^V ^p" = 0, which 
is the same as the description in the standard Einstein 
gravity, where the rest inertial mass of particle is con- 
served. 

Any change in the gravitational mass in the Jordan 
frame inevitably changes the gravitational field, which, in 
turn, will change the inertial frame. Any change in the in- 
ertial frame will induce the well-known "frame-dragging" 
effect which was recently confirmed by the Gravity Probe 
B mission 34]. The inertial "frame-dragging" effect 
can be equivalently considered as the "inertial mass- 
dragging" effect by assuming that the inertial frame is 
unchanged while the inertial mass of particles defined in 
this frame is rescaled. This understanding can help us 
interpret that after the conformal transformation, in the 
Einstein frame, the equation of motion for particles in 
the inertial frame is governed by the varying mass equa- 
tion of motion u^V ^p v = md** In yJ~F = ^(jj)", which 
is one of the well established physics in General Relativ- 
ity [26J, 35]. The rest inertial mass in the Einstein frame 
is indeed rescaled as fhi = m-jm an d is no longer con- 
served. 

Considering the Mach principle, there should be no 
difference among gravitational theories described in dif- 
ferent frames. Thus a viable f(R) model in the Jordan 
frame should be consistent with the viable model in the 
Einstein frame in describing cosmology. The result ob- 
tained in the work has solid physics root. 
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